STABILITY OF PID-CONTROLLED LINEAR TIME-DELAY 
FEEDBACK SYSTEMS 
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Abstract. The stability of feedback systems consisting of linear time-delay 
plants and PID controllers has been investigated for many years by means 
of several methods, of which the Nyquist criterion, a generalization of the 
Hermite-Biehler Theorem, and the root location method are well known. The 
main purpose of these researches is to determine the range of controller param- 
eters that allow stability. Explicit and complete expressions of the boundaries 
of these regions and computation procedures with a finite number of steps 
are now available only for first-order plants, provided with one time delay. In 
this note, the same results, based on Pontryagin's studies, are presented for 
arbitrary-order plants. 



1. Introduction 

The feedback structure considered in this note is depicted in Fig. 1 and the 
related transfer functions of the process P{s) and the controller C(s) are given by 



(1.1) 



(1.2) 



P(s) 



Pais) 



C{s) = i^p + — + KdS, 



where K is the plant steady-state gain, Ti and Zi the plant time constants, L is 
the positive plant time delay and Kp, Ki and Kd are the parameters of the PID 
controller. 
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Figure 1 . Feedback control system 



Complete explicit expressions of the boundaries of the stability regions in first- 
order plants have been found in [1] with a version of the Hermite-Biehler Theorem 
derived by Pontryagin, in ^2^ with the Nyquist criterion, and in [3] with the root 
location method. Moreover the second-order plants have been investigated in ^ by 
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means of a graphical approach; the results obtained are correct, but the stability 
conditions are not all explicit and no finite number of required computation steps 
is specified. Finally arbitrary-order plants have been studied with the Nyquist 
criterion in [5] , but P and PID controllers with a given Kp separately are considered 
and no information about the set of the process parameters that allow stability is 
given. 

This note can be considered as an extension of jlj to the arbitrary-order plants 
and is organized as follows. In Section 2 all the analytical expressions that will be 
used in the next sections are evaluated in detail. In Section 3 a process transfer 
function without zeros is considered and the stability regions are explicitly evaluated 
by means of a version of the Hermite-Biehler Theorem derived by Pontryagin, 
already used in [I]. A second-order plant is studied as example and the related 
stability regions are determined and plotted in two figures. In Section 4 a process 
transfer function with zeros is considered and the stability regions are found by 
means of a new theorem. The two procedures of the Sections 3 and 4 are essentially 
equal, consist of a finite number of steps and yield the stability regions in both 
process and controller parameters planes. In Section 5 some conclusive remarks are 
given. 

The importance of explicit expressions of the boundaries of the stability zones 
has been enhanced by the introduction of the controller tuning charts in [6J (used 
also in ff]). 



2. Preliminaries 
The closed-loop transfer function T(s) of the system is given by 

^ ' ' ^ ' sPd{s)eL' + K{K, + Kps + KdS^)Pn{s) ' 

According to the Pontryagin's studies, presented in ^ and summarized in [9], it is 
necessary that T{s) has a bounded number of poles with arbitrary large positive 
real part for stability. This holds if the denominator of T{s) has a principal term 
apqS^e"^" (in our case, where p = n -\- 1 and g = 1, it exists if m <— n — 1) and the 
function Xpi^)j coefficient of s^, (in our case Xp(s) = ^'^^YilZi'^i for 77i < rt — 1 
and Xp(s) = e'' ^ 11^=1 Ti + K Kd ni=i" f^^' ra ~ n — 1) has all the zeros in the 
open left half plane. This happens if one of the following conditions is satisfied: 

(a) m < n — 1 

(b) m^n-ltind\KKd 11 -l'" z}<\ n^l" . 

The denominator of T{s), given by (|2.ip . divided by P„(s)/L and hence named 
H{s), can be written, according to (|l.ip . as 

(2.2) i/(s) = L.e^^^il±M +LKiK, + KpS + Kds'). 

n,^i (1 + z,s} 

Since all the poles of the closed-loop transfer function T(s) are zeros of H{s) and 
a system is stable if no pole of T{s) lies in the right half-plane, the above system is 
stable if no zero of H{s) lies in the right half-plane. For process transfer functions 
without zeros, examined in Section 3, H(s) is a quasi-polynomial and a version 
of the Hermite-Biehler Theorem derived by Pontryagin is employed. For process 
transfer functions with zeros, examined in Section 4, H{s) is a quasi-polynomial 
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divided by a polynomial and a new theorem, proved by use of the Principle of the 
Argument, is employed. 

Now, before the explanation of the proposed procedures, let us evaluate all the 
expressions that will be used in the next sections. It is convenient to introduce the 
normalized time referred to the plant time delay L and the dimensionless parameters 
a = L s, ti = Ti/L, Zi = Zi/L, h — K Kp, hi = K KiL and hd = K Kd/L, in order 
to obtain equations independent of the real values of the parameters. Applying 
these simplifications, (|2.2p becomes 



(2.3) 



Hia) 



ae 



.n::i(i+t.-) 



hi + ha + hdcr^ 



n::r(i+-«'^) 

Moreover, assuming Pd{j y/L) = A{y) +jB{y) and Pn{j y/L) = C{y) + jD{y), the 
real and the imaginary components F{y) and G{y) of H{a), calculated for a = jy, 
are given by 

F{y)^K-F^{y) 



(2.4) 

(2.5) 

where 

(2.6) 

(2.7) 

(2.8) 



(2.9) 



(2.10) 



(2.11) 



G{y)^y[h-G^{y)] 

Fi{y) = y[Q{y)cos{y) + P{y)sin{y)] 

Gi{y) = -P{y)cos{y) + Q{y)sin{y) 

he^hi- hd 
A{y)C{y) + B{y)D{y) 



Piy) = 
Q{y) 



C^{y) + D^y) 
-A{y)D{y) + B{y)C{y) 
CHy) + D^y) 



i—int{n/2) 

A{y) = l+ J2 U{n,2t){-l)V' 
1=1 

i=mi((n-l)/2) 

i—int{m/2) 

C{y) = l+ n"^,2*)(-l)V' 



(2.12) 



Diy) 



i=int((m-l)/2) 

Y n™,2i+i)(-i)V'+i. 



1=0 



For sake of clarity, U and V are the symmetric expressions of the time constants ti 
and Zi; U{n, k) is the sum of the products of k different U selected among the 
total n (for example 17(3, 2) = tit2 + ^2^3 + ^3^1 and [7(3, 3) = ti<2*3)- 
The derivative of G{y) with respect to y is given by 



(2.13) 



G'{y)^h-G,{y)-yG\{y). 
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Assuming G{y) — and h = —1 in (j2.5p . one obtains tan{y/2) — E{y) where 



(2.14) Eiy) . ~Qiy)^VP'(y) + Q'(y) 



(2-17) 



1 + Piy) 

It is easy to check that E{y) exists for y — only if X)i=i ^1 > Si^™ -^i^- The 
derivative of E{y) with respect to y, evaluated a,t y — and named -E'(O), is given 

by 

E'{0) = +0.5(-L/(n, 1) + T/(to, 1)) 

(2.15) , 

± 0.5 1) - 2 U{n, 2) - ^^(m, 1) + 2 V{m, 2). 

Denoting by E_{y) and E^{y) the two branches of £'(?/) related respectively to 
the minus and plus signs, their derivatives E'_{0) and E'_^_{0) are higher than the 
derivative of tan{y/2), equal to 0.5, depending on <I>i and $2, given by 

(2.16) $1 = l + f/(n, 1) -y(m, 1) 
$2 = + 1 + 2 U{n, 1) + 2 U{n, 2) - 2 U{n, l)V{m, 1) 

+ 2 V^{m, 1) - 2 V^(m, 1) - 2 y(m, 2). 

In detail, the number of the derivatives E'_{0) and E'^{0) higher than 0.5 are the 
following: zero if $1 > and <i>2 > 0, one if $2 < 0, and two if $1 < and $2 > 0. 
Let us denote by Ed 

(2.18) Ea = taniyt/2) - E{yt) 

where yt, corresponding to equal derivatives with respect to y of tan{y/2) and E{y), 
is a root of 0.5(1 + tan^{yt/2)) = E'{yt). 

Differentiating (|2.7p with respect to y once and twice, one obtains 

(2.19) = + Q{y))cos{y) + (P(y) + Q'{y))sin{y) 
ay 

(2 20) ^ + ^"^^^ + ^ Q\y))cos{y) 

+ {~Q{y) + Q"{y) + 2P'(y))sin(y). 

It is worthwhile to note that Gi(0) = -1, dGi{0)/dy = 0, and (PGi{0)/dy^ = $2, 
where $2 is given by (|2.17p . Evaluating hm = \Gi{ym)\, where is a root of 
dGi (y) / dy given by (|2.19p , one obtains 

.2 21) _ |^'(2/m) + Q'(y„0 + Piym)Q'{ym) ~ P'iym)Q{ym)\ 

Eliminating cos{y) and sinly) from F{y) = and G(y) = 0, given by (|2.4p and 
5]), yields 



(2.22) /le = y^P^W+QHjj)^. 

Denote by Tai and Ffci the two straight lines whose equations in the (/i^, /i(i)-plane 
are obtained introducing in (|2.8p respectively y = ^ = ^eoi and y = y^i, 
^ — hebi (see Figs. 2 and 5); denote further by K:, Ui and the vertices of a 
triangle, whose sides are the axis hd and the two lines Tai and Tbi. The coordinates 
of these vertices are given by 

(2.23) h,{Vi) = ~ ^1^'^'" ; /id(KO = _+ 

2/f)i ^ai ^fci 2/ai 
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(2.24) ham - ^heu/yl; h4W,) = ~hea^/yl,■ 

Considering (I2.22p . the coordinates hd{Ri) and hd{Si) of the points lying on Fq 
and Tbi at hi = hi{Vi) are given by 



(2.25) /id(-Ki) = 2 — - sign[hebi\ 

Vbi Vbt 



(2.26) /id(5i) = — 2 — ~ sign[heai\- • 

Vai Vo-i 

It is easy to check that, when i = oo, the absolute values of hd{Ui) and hd{Wi) , if 
rn < n — 1, are equal to oo and, if m = n — 1, to Hd{oo) given by 

(2.27) Hd{oo) ^\U{n,n)/V{m,m)\. 

3. Process transfer function without zeros 

When the process transfer function does not have zeros, m — and thus C{y) = 
1, D{y) = 0, P{y) = A{y), Q{y) — B{y) hold; therefore, the function H{a), given by 
(|2.3p . is a quasi-polynomial and the Pontryagin's results are integrally applicable. 
The following two conditions derived from Theorem 3.2 of jTj and from Theorem 
13.7 of [9j, respectively, must be satisfied in order to have a stable system: 

• Condition no. 1 

Consider that the principal term of H{(j), given by (|2.3p . is o'"+^e°', set 
H{j v) = F{y) + jG{y) and let e be an appropriate constant such that the 
coefficient of y"^^ in G{y) does not vanish at y = e. The number Nr of the 
real roots of G{y) in the interval —2rTT + e<y<2rTT + e for sufficiently 
large r must be 

(3.1) Nr = Ar + n + l. 

• Condition no. 2 

For all the zeros y = yo of the function G{y) the inequality G' {yo)F{yo) — 
G{yo)F'{yo) > 0, that is G'{yo)F{yo) > 0, must hold. 
In order to study both stable and unstable free-delay plants, the following two 
cases, adopted also in [T], are considered: 

(1) U{n, n) > (even number of negative plant time constants Ti) 
h > —1 and hi > 0. 

(2) U (n, n) < (odd number of negative plant time constants Ti) 
h < —1 and hi < 0. 

From (|2.5p . (|2.9p and (|2.10p it follows that the coefficient of the highest degree 
of y in G{y) is U{n,n)cos{y) when n is even and U{n,n)sin{y) when odd; hence 
we assume e = if n is even and e = O.Stt if odd. 

Typical functions Fi{y) and Gi(y) are plotted in Fig. 2; according to p.Sp 
there is one root of G{y) at y = and one for each intersection of Gi(y) with the 
horizontal line having the ordinate equal to a given h. 

Denoting by the number of the intersections between Gi{y) and h = —1 
corresponding to N^. in Fig. 2 and assuming that no local minimum or maximum 
of Gi{y) is equal to —1 for y 7^ 0, the relationship between and is given by 

N^^Nr-l = Ar + n for C/(n,n)$2<0 

Ne^Nr~3 = Ar + n-2 for C/(n,n)$2>0, 




Figure 2. Typical plots of and Gi{y) 



where $2 is according to (|2.17|) : (|3.2|) can be easily checked considering that from 
(lOl . ({239)) and (12:201) it follows Gi(0) = -1, G[iO) = and G'{{0) = *2, and 
also that h > —1 must hold if U{n,n) > and /i < — 1 if U{n,n) < 0. Since 
/i = — 1 is the common limit value of h for the two considered cases {U{n,n) < 
and U{n,n) > 0), the existence of the intersections given by (|3.2p represents a 
prerequisite of a plant to be made stable. 

The number can be evaluated by counting the intersections of the plots of 
tan{y/2) and E{y), given by (|2.14p . From (|2.14p it follows that E{y) is an odd 
function of y; moreover, assuming 5*^ = sign[A{+co)] and Sb = sign[B{+co)], 
i?_(+oo) and i?+(+cx)) can be expressed as 

• n even 

E^{+oo) = -Sa 1 and E+{+oo) = +S'a 1 . 

• n odd 

S_(+oo) -~Sa 00 for S'h > and £;_(+oo) = -Sa for Sb < 0, 
i;+(+oo) = +Sa for > and E+{+oo) = +S'a 00 for Sb < 0. 
If £'(?/) has no pole, splitting into A^ei (|y| < tt) and Ne2 and considering the 
above described behavior of E{y) at y = 00 and also at y = 0, one obtains 

• |y| < 71" 

iV^i = if $1 > and $2 > 0, 

iVei = 2 if $2 < 0, 

iVei = 4 if $1 < and $2 > 0, 

where $1 and $2 are given by (|2.16p and (|2.17p . 

• —2 rn + e < y < — tt; 7r<y<2r7r + e (see Fig. 3 (a)) 
Ne2 = 4r — 2 if ri is even, 

iVg2 = 4r - 1 if n is odd and A{+oo)B{+oo) > 0, 
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Ne2 = 4r - 3 if 71 is odd and A[+oo)B{+oo) < 0, 

where sign[A{+oo)B{+(X))] = sign[—{—l)'"-U{n,n~ l)U{n,n)] . 

It is clear that, if E{y) has no pole, Ne will be always lower than the value required 
by (|3.2[) for enough large n. A positive solution can be reached only if E{y) is 
provided with a suitable number of poles, since A^e2 is increased by one for each 
added pole (see Fig. 3); this happens if -E^ > for case (bl), if Ed < for case (b2) 
and without further condition for case (b3), where Ed is given by ()2.18p . Since the 
denominator of E{y) is a polynomial of y'^ of degree d = int{n/2), the maximum 
number of poles of E{y) is equal to 2 d and the actual number can be determined 
by means of the Sturm Theorem, as detailed in Appendix A. 




Figure 3. Plots of tan{y/2) and E{y) for y > tt 



The sought-after procedure can be detailed as follows: 

(1) Process parameters (see Fig. 4) 

The stability region in the process parameters plane is that where p.2p 
holds; its boundary line is a proper set of the boundary lines of the zones 
with different numbers A''e, that is, <I>i = and $2 = according to (|2.16[) 
and (|2.17p and "ifi j = 0, as explained in Appendix A, and eventually Ed = 
given by (|2.18p . Since the expressions of these boundary lines are functions 
of ti = Ti/L, it is possible to evaluate the stability range of L when the 
parameters Ti are known and, conversely, of each Ti when L and the re- 
maining Ti are known. Moreover, if one needs only to know if a given plant 
can be made stable, it is not necessary to determine the stability regions, 
but it is enough to examine the plots of E{y) and tan{y/2) and to check 
whether the number of the intersections Ne satisfies p.2|) . 
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(2) Controller parameter h (see Fig. 2) 

The requirement, stated as Condition no. 1, is fulfilled if the selected value 
of h is included in the interval from —1 to hp for U{n, n) > or from /i„ to 
— 1 for U{n, n) < 0, where hp and hn are respectively the local maxima and 
minima of Gi{y) nearest to —1; both can be calculated by introducing in 
(|2.7p the related root of dGi{y)/dy obtained from (|2.19p . A finite number 
of these maxima or minima must be examined in order to find the nearest 
to —1, exactly up to the first value of j/„i higher than y^i- This limit y^i is 
the largest positive root of the derivative with respect to ym of hrm since 
h„n given by (j2.2ip . monotonically increases for ym > yri- It is obvious 
that, if such root does not exist, only the first value must be considered. 

(3) Controller parameters hi, hd (see Fig. 5) 

Considering (|2.4|) . (j2.8|) and (|2.13|) . the requirement, stated as Condition 
no. 2, is fulfilled if the following inequalities 

hi - hdvl < Fi{yo) G;(yo) > 

K-hayl>F^{ya) if G'M <0 

hold for each root y = yo oi G{y), given by (|2.5p . evaluated with a value of 
h included in the above specified interval. 

The stability region in the {hi,hd)-p\ane consists of the intersection of 
a finite number of triangles; each of them is related to a couple of roots 
yai and y^ of G{y), has the axis hd and the two straight lines given by 
p.3p as sides and the points Ui, Vi and Wi as vertices, whose coordinates 
are given by (|2.23p and (|2.24p (see Figs. 2 and 5). Since, as i — > oo, 
hd{Ui) and hd(Ri) approach ±oo and hdiyVi) and hd{Si) approach Too, 
each triangle includes definitely the first one when yu > yr2- This limit yr2 
is the bigger among y^i and the largest root of the derivatives with respect 
to yu of hd{Ui) and hd{Ri), given by (|2.24p and (|2.25p . Therefore, it is not 
necessary to examine the triangles for yu > j/,-2- 
A second-order plant, whose transfer function is without zeros, is considered 
as example and the stability regions of the plant and the controller parameters are 
depicted respectively in Figs. 4 and 5. In Fig. 4 this region consists of the following: 

• Zi: U{n, n) > 0; $1 > and $2 > 

• Z2: u\n, n) > 0; $1 < and $2 > 

• Z3: U{n,n) < 0; $2 < 

where $1 1 + ii + <2 and $2 = 1 + 2fi + 2t2 + 2tit2 according to and 
(|2.17p . The required number of the intersections between E{y) and tan{y/2), 
equal to 4 r as per (|3.2p . coincide with the actual number only in this zone. 

Since P{y) = 1 — ^1^22/^ Qiu) = (^i + ^2)2/, considering the point (ti ~ 0.6; 
t2 — 0.8) lying in Zi of Fig. 4, one obtains yp — 1.778 from (|2.19p for the first root 
y = yp oi dGi{0)/dy and hp = Gi{yp) = 2.330 from (|2.7p . Since —l<h<hp must 
hold, let us assume h equal to 0.5; for the first two roots higher than zero of G{y) 
one obtains yai = 0.863 and yhi = 2.498 from (|2.5p and hence /leai — ^i(2/ai) = 
1.099 and hebi = Fi(y6i) = -9.985 from dSl]). From ^jT^ and f?:^ it follows 
h,{Vi) = 2.600, hdiVi) = 2.016, hd{Ui) = 1.600 and hdiWi) = -1.476. Similarly, 
for the third and fourth roots of G{y) one obtains ya2 = 5.285 and yb2 — 8.191, 
Ka2 = 76.290, heb2 = "272.288, hd{U2) = 4.058 and hdiW2) = -2.732; moreover 
hdiR2) = 4.097 and /id(S'2) = -2.638 from fiTmi and (^7^ . 











Zi 


r 


1 \ 










Z2 \ 




\ Z3 



h 

Ti : 1 + + t2 =0 

r2 : l + 2(ti + ^2) + 2tii2 = 

Figure 4. Stability zone of the process parameters for n = 2 and 

Pn{s) = 1 




Figure 5. Stability zone of the controller parameters hi and hd 
for 71 = 2, P„(s) = 1, 0.6, t2 = 0.8, h ^ 0.5 



4. Process transfer function with zeros 
The function H{a), given by (|2.3p . can be rewritten as 

(4.1) ^(-) = -""^^"^-^^^^' 

v (to, to) 1 + Hd 

where 

^^U(n,n — i) _n , _i , , F(m,n— 1 — z) 



X ^ 1/(to, m — i) 
V{m,m)(T^ 



m<n— 1 or m^n— 1 and |/idT^(TO, m)| < |C/(n,n)|, as explained in Section 2. 
Since the denominator of H{a) is the numerator of the plant transfer function, the 
poles of H{a) are the zeros of the plant transfer function {ai ~ —1/zi). In this case 
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the proposed procedure will be in accordance with Theorem 4.1, a generalization 
of the theorems applied in Section 3 and of Theorem 13.5 of 9J; it will be here 
enunciated and proved by use of the Principle of the Argument. 

Theorem 4.1. Consider a function H (a) of the form in setH(jy) — F(y) + 

jG{y) and let e be an appropriate constant such that the coefficient of in 
the numerator of G{y) does not vanish at y — e. Assume further that the rUp poles 
with positive real part of H(a) lie all in the rectangle R, described by the inequalities 
— 2r7r + e <= y <= 2r7r + e, a;>0 (a = x + j y). Suppose finally that the function 
H{a) does not assume the value of zero on the imaginary axis, that is, H{jy) ^ 0. 
All the zeros of H{(7) lie to the left of the imaginary axis if and only if: 

(a) The vector w = H{j y) for real y ranging from —oo to +oo continually 
revolves in the positive direction at a positive velocity, that is, the inequality 
G"(yo)^(yo) < is satisfied for each root y = yo of G{y). 

(b) The number Nr of the roots of G{y) in the interval —2rii + e y <= 
2r7r + e for sufficiently large r is 

(4.2) iVr = 4r + n + 1 - m + 2mp. 

Proof. Denote by Ci , C2 , C3 and C4 the corners of the rectangle R whose coordi- 
nates in the cr-plane are respectively cri = + 2r7r + e), (T2 = -l-(X)-|-j(— 2r7r-|-e), 
0-3 — +00 + j{+2r'K + e) and 0-4 = + j{+2r'K + e); the arguments 0{ai) and 
0{a4) of H{a) are given by 0{ai) = — 2r7r + e — 0.5(ri + 1 — m)7r + rj + Si and 
0((T4) = +2r7r + e + 0.5(n + 1 — m)7r -\- rj + 64, where ry = if V{m, m)/U{n, n) > 
or ?7 = TT if V{m,m)/U{n,n) < 0, and 5i — > and 64 ^ simultaneously with 
1/r. Denote by the number of the zeros of H{a) lying in the rectangle R, by Va 
the variation of the argument of H(a) in the counterclockwise direction as a moves 
around the contour of R from Ci to C4 through C2 and C3, and by Vf, as a moves 
directly from C4 to Ci. Using the Principle of the Argument yields 

(4.3) Va + Vb = {N,-mp)2 7:. 

Since A^^ = for stability and , for r — > 00, Vf, = —Nrir as per condition (a) and 
Va = 6*4 - 6*1 = +4r7r + (n + 1 - m)n, ([42]) follows from ([43|) and the condition 
(b) is satisfied. □ 

The procedure detailed in Section 3 for process transfer functions without zeros 
is fully applicable to process transfer functions with zeros; it is only necessary 
to replace ()3.ip with (|4.2|) and to consider C{y) and D{y) as functions given by 
p. lip and (|2.12p instead of C{y) = 1 and D{y) = 0. Moreover, for m = n — 
1, the rectangle, according to (|3.3p for j/o — +00 and provided with horizontal 
sides symmetric with respect to the axis hi at a distance given by (|2.27p . must be 
considered in the (hi, ft,(j)-plane. 

5. Conclusions 

In this note, both stable and unstable delay-free arbitrary-order plants, provided 
with one time delay and PID controller, have been examined and the related stabil- 
ity regions in process and controller parameter spaces have been determined by use 
of the Pontryagin's studies. The proposed procedure, consisting of a finite number 
of steps, yields explicit expressions of the boundaries of the stability zone for the 
controller parameters. These results can be implemented in tuning charts, which 
become a complete tool for the design and the maintenance of control systems. 



11 



Appendix A. Sturm Theorem 

The Sturm Theorem states that the number of real roots of an algebraic equation 
with real coefficients whose real roots are simple over an interval, the endpoints of 
which are not roots, is equal to the difference between the numbers of sign changes 
of the Sturm chains formed for the interval ends. 

Given a function f{x) = fo{x) of degree d, assume fi{x) = df(){x)/dx and define 
the Sturm functions by 

Mx)^-h-2{x)+.h-l{x)\^^P^ , 

where [fi-2{x)/fi-i{x] is a polynomial quotient. These functions can be written as 

j=d-mt((i+l)/2) 

f,{x)= ^^r^' 0<i<2d-l, 

j=o 

where ipij depends on the coefficients of x in f{x). 

In our case a; = and/o(a;) = C^{^/x)+D^{^)+A{^)C{y^)+B{^)D{^/^) 
hold; since the roots must be positive, the required interval is from x = to x = +00 
and, therefore, the signs of each Sturm function at these ends are the signs of ipij 
evaluated respectively for j = and j ^ d~ int{{i + l)/2). 
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